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Abstract 



The recent progress in establishing the existence of finite neutrino masses and mix- 
ing between generations of neutrinos has been remarkable, if not astounding. The 
\ combined results from studies of atmospheric neutrinos, solar neutrinos, and reactor 

I antineutrinos paint an intriguing picture for theorists and provide clear motivation 

' for future experimental studies. In this review, we summarize the status of experi- 

, mental and theoretical work in this field and explore the future opportunities that 

\ emerge in light of recent discoveries. 
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1 Introduction and Historical Perpective 



For almost 70 years, neutrinos have played a pivotal role in the the quest 
to understand elementary particles and their interactions. The neutrino was 
actually the first particle proposed by a theorist to guarantee the validity 
of a symmetry principle: Pauli boldly suggested [1] that neutrinos (invisible 
then and remaining so for another twenty five years) are emitted together with 
electrons in nuclear beta decay to salvage both energy and angular momentum 
conservation in the beta decay process. Experimental confirmation of Pauli's 
hypothesis required many decades of experimental effort, but the discovery 
of the antineutrino in the fifties by Reines and Cowan [2] and subsequent 
experiments in the early sixties ultimately led to the Nobel Prizes for Reines 
(1995) and Lederman, Schwartz and Steinberger (1988). More recently, the 
2002 Nobel prize was awarded to Davis and Koshiba for their seminal roles 
in the development of neutrino astrophysics (along with Giaconni for x-ray 
astrophysics) . 
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The standard model of electroweak interactions, developed in the late 1960's, 
provided a theoretical framework to incorporate the neutrinos as left-handed 
partners to the charged leptons, organized in two generations along with the 
quarks. The subsequent discovery of charmed quarks and the third generation 
of quarks and leptons completed the modern view of the standard model of 
electroweak interactions. This version possessed additional richness to incor- 
porate CP violation, and further efforts to unify the strong interaction led to 
the development of Grand Unified Theories (GUTs). These GUTs provided a 
natural framework for nucleon decay and for neutrino masses, and motivated 
many experiments in the field. 

Following on the successes of big-bang nucleosynthesis (BBN) and the discov- 
ery of the cosmic microwave background, it also became clear that neutrinos 
were potentially major players in the history of the early universe. The in- 
creasing evidence for the existence (and then dominance) of dark matter in 
the universe then led to the economical and seductive hypothesis that neutri- 
nos, with small but finite mass, could provide the mass to explain the dark 
matter. This set the stage for a major experimental assault on the issue of 
neutrino mass and its role in cosmology, and provided substantial impetus to 
a worldwide program of experiments addressing the issues of finite neutrino 
mass and the possibility of mixing between generations. 

Although it now appears that neutrinos are not likely the source of dark matter 
in the universe, the experimental evidence obtained in the last decade for finite 
neutrino masses and mixing between generations is strong and irrefutable. The 
pattern of masses and mixing angles emerging from the experiments provides 
an intriguing glimpse into the fundamental source of particle mass and the role 
of fiavor in the scheme of particles and their interactions. The scale of neutrino 
mass differences motivates new experimental searches for double beta decay 
and end-point anomalies in beta decay, as well as new studies of oscillation 
phenomena using accelerators, nuclear reactors, and astrophysical sources of 
neutrinos. 

In this review, we attempt to synthesize these themes, present a concise and 
coherent view of the experimental and theoretical developments leading to the 
current picture, and motivate the future explorations necessary to resolve the 
remaining issues ^ . We begin with a summary of the theoretical motivation 
for studying neutrino masses and mixing, along with a development of the 
phenomenological framework for interpreting the experiments. This includes 
neutrino oscillations in vacuum as well as in matter, CP violation, ordinary 
beta decay, and double beta decay. We then review the recent experimental 



The reference list is, for reasons of brevity, far from complete. We apologize to 
those whose work is not quoted, and refer to numerous reviews and monographs to 
provide a more complete reference list. 
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data that contribute to our present knowledge of these neutrino properties. 
We discuss the successful neutrino oscillation measurements, including the key 
contributions from various solar neutrino experiments, cosmic ray-induced at- 
mospheric neutrino studies, and the recent dramatic results from the Sudbury 
Neutrino Observatory and the KamLAND reactor neutrino experiment. Im- 
portant additional information is obtained from the experimental attempts 
that have thus far only yielded limits on neutrino masses such as double beta 
decay and tritium beta decay. We also will briefly discuss the role of massive 
neutrinos in cosmology, and the corresponding constraints. 

The discussion of the future experimental program is separated into two dif- 
ferent time scales. In the near-term future, the planned experiments will focus 
on a better determination of the mixing matrix parameters, in particular ^13, 
and resolution of the LSND puzzle (right or wrong). Further studies in the 
longer-term will of course depend on the outcome of these measurements, but 
there is substantial interest in pursuing the search for CP violation, the pos- 
sibility that neutrinos are Majorana particles and other issues related to the 
mechanism(s) responsible for the observed phenomena. 



2 Theoretical Framework 

It is natural to expect that neutrinos have nonvanishing masses, since they are 
not required to be massless by gauge invariance or other symmetry principles. 
Moreover, all other known fermions, quarks and charged leptons, are massive. 
Nevertheless, neutrinos are very light, much lighter than the other fermions, 
and this striking qualitative feature needs to be understood, even though we 
do not know why e.g. the electron mass is what it is known to be, and why 
muons and taus are heavier than electrons. 

The most popular explanation of the small neutrino masses is the "see-saw 
mechanism" [3] in which the neutrino masses are inversely proportional to 
some large mass scale M^. Simply put, if neutrinos are Majorana particles, 
i.e., indistiguishable from their own antiparticles, they have only two states, 
corresponding to the two possible spin orientations. On the other hand, Dirac 
particles, distinct from their antiparticles, have four states, two spin orienta- 
tions for the particle and antiparticle. In the "sec-saw mechanism" there are 
still four states for each neutrino family, but they are widely split by the "Ma- 
jorana mass term" (explained later) into the two-component light neutrinos 
with masses and the very heavy (sterile) two-component neutrinos with 
the mass M->^, in such a way that McM-^ = M|,. Further, if we assume that 
the "Dirac mass" Mp is of the same order of magnitude as the Dirac fermion 
masses (masses of quarks and charged leptons), we can understand why Mc is 
so small, provided M-^ is very large. If that explanation of the small neutrino 
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mass is true, then the experimentally observed neutrinos are Majorana parti- 
cles, and hence the total lepton number is not conserved. Observation of the 
violation of the total lepton number conservation (we explain further why it is 
difficult to observe it) would be a signal that neutrinos are indeed Majorana 
particles. 

Neutrinos interact with other particles only by weak interactions (at least as 
far as we know); they do not have electromagnetic or color charges. That 
further distinguishes them from the other fermions which also interact elec- 
tromagnetically (charged leptons and quarks) and strongly (quarks). 

It is well known that the weak charged currents of quarks do not couple to a 
definite flavor state (or mass state), but to linear combinations of quark states. 
This phenomenon is described by the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix, which by convention describes linear combinations of the charge — e/3 
quark mass eigenstates rf, s, b. The mixed states d', s', b', form weak doublets 
with the u, c, t quark mass eigenstates. Consequently, the three generations 
of quark doublets are not independently preserved in charged-current weak 
processes. 

An analogous situation is encountered with neutrinos. There, however, one 
can directly observe only the 'weak eigenstates', i.e., neutrinos forming weak 
charged currents with electrons, muons or tau ^ . Information on the neutrino 
mixing matrix is best obtained by the study of neutrino oscillation, a quantum 
mechanical interference effect resulting from the mixing. As a consequence, 
like in quarks, the flavor (or family) is not conserved and, for example, a 
beam of electron neutrinos could, after travelling a certain distance, acquire 
a component of muon or tau flavor neutrinos. We describe the physics and 
formalism of the oscillations next. The formalism has been covered in great 
detail in numerous books and reviews. Thus, we restrict ourselves only to 
the most essential points. The interested reader can find more details in the 
monographs [4,5,6,7,8,9,10] and recent reviews [11,12,13,14,15,16,17,18]. 

2.1 Neutrinos in the standard electroweak model 

In the standard model individual lepton charges (Le = 1 for e" and and 
Lg = —1 for e"*", and Pg and analogously for and L^) are conserved. Thus, 
processes such as //"^ — > e"*" + 7, or Kl + //^ are forbidden. Indeed, 

^ The analogy is not perfect. For quarks, the 'mass eigenstates' are also 'flavor 

eigenstates', with labels d, s, b. Since the CKM matrix is nearly diagonal, the labels 
d' , s' , b' can be used for the 'weak eigenstates'. In neutrinos, on the other hand, only 
the 'weak eigenstates' have special names, v^^^h^^t while the mass eigestates are 
labelled I'l, 1^2, 1^3 and cannot be associated with any particular lepton flavor. 
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such processes have not been observed so far, and small upper limits for their 
branching ratios have been established. 



Based on these empirical facts, the standard model places the left-handed 
components of the charged lepton and neutrino fields into the doublets of the 



while the right-handed components of the charged lepton fields are singlets. 
The right-handed components of the neutrino fields are absent in the standard 
electroweak model by definition. 

As a consequence of this assignment, neutrinos are deemed to be massless, and 
the individual lepton numbers, as well as the total one, are strictly conserved. 
Note that with these assignments neutrino masses do not arise even from loop 
corrections. Thus, observation of neutrino oscillations, leading to the neutrino 
flavor nonconservation, signals deviations from this simple picture, and to the 
'physics beyond the standard model'. 

Note also that studies of e+e~ annihilation at the Z-resonance peak have 
determined the invisible width of the Z boson, caused by its decay into un- 
observable channels. Interpreting this width as a measure of the number of 
neutrino active flavors, one obtains Ni, = 2.984 ± 0.008 from the four LEP 
experiments [20]. We can, therefore, quite confidently conclude that there are 
just three active neutrinos with masses of less than Mz/2. (The relation of this 
finding to the fact that there are also three flavors of quarks is suggestive, but 
so far not really understood.) Besides these three active neutrino flavors there 
could be other neutrinos which do not partipate in weak interactions. Such 
neutrinos are called 'sterile'. In general, the active and sterile neutrinos can 
mix, and thus the sterile neutrino can interact, albeit with a reduced strength. 

Big Bang Nucleosynthesis (BBN) is sensitive to the number of neutrino fla- 
vors which are ultrarelativistic at the 'freezeout' when the n <-> p reactions 
are no longer in equilibrium and therefore it is perhaps sensitive to (almost) 
sterile neutrinos. However, analysis of BBN (i.e., of the abundances of *^He, d, 
^He, and ^Li) are consistent with 3 neutrinos, and disfavors A^,^ = 4 for fully 
thermalized neutrinos [19]. 



group SU{2)l, 
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2.2 Neutrino mass terms 

In a field theory of neutrinos the mass is determined by the mass term in the 
Langrangian. Since the right-handed neutrinos are absent in the standard elec- 
troweak model, one can either generalize the model and define the mass term 
by using the ideas of the see-saw mechanism, or one can add more possibilities 
by adding to the three known neutrino fields new fields, corresponding to 
possibly heavy right-handed neutrinos i/jR. The mass term is then constructed 
out of the fields ly^L, their charge-conjugated and thus right-handed fields (f^i)"^ 
and from the fields VjR, {i^jRY- There are, in general, two types of mass terms. 



(The term with y'^e,iL^e,L is left out assuming that the corresponding coefficients 
are vanishing or negligibly small.) 

The first term is a Dirac mass term, analogous to the mass term of charged 
leptons. It conserves the total lepton number, but might violate the individual 
lepton flavor numbers. 

The second term is a Majorana mass term which breaks the total lepton 
number conservation by two units. It is allowed only if the neutrinos have no 
additive conserved charges of any kind. 

If there are 3 left handed neutrinos vil, and n additional sterile neutrinos VjR 
it is convenient to define 



The matrix is a symmetric complex matrix and M.^ is the transposed 
matrix. After diagonalization it has 3 -h n mass eigenstates Vk that represent 
Majorana neutrinos {v^ — Uk). 

From the point of view of phenomenology, there are several cases to be dis- 
cussed. The seesaw mechanism corresponds to the case when the scale of M.'^ 
is very large. There are three light active Majorana neutrinos in that case. 

On the other hand, if the scale of Ai^ is not too high when compared to the 
clcctroweak symmetry breaking scale, there could be more than three light 
Majorana neutrinos, mixtures of active and sterile. Finally, when A4'^ = 
there are six massive Majorana neutrinos that merge to form three massive 




(2) 
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Dirac neutrinos. The unitary matrix diagonalizing the mass term is a 3 x 3 
matrix in that case. In these latter cases the 'natural' explanation of the 
lightness of neutrinos is missing. 

We shall not speculate further on the pattern of the mass matrix, even though 
a vast literature exists on that subject (see e.g. [13] for a partial list of recent 
references). In the interesting case of light left-handed Majorana neutrinos, 
the mass matrix is determined by neutrino masses, N{N — l)/2 mixing 
angles, [N — 1){N — 2)/2 CP violating phases common to Dirac and Majorana 
neutrinos, and {N — 1) Majorana phases that affect only processes which 
violate the total lepton number (for the discussion of these Majorana phases 
see e.g. [21]). 

2.3 Neutrino oscillation in vacuum 

As stated earlier, the neutrinos participating in the charged current weak in- 
teractions (the usual way neutrinos are observed and the only way their flavor 
can be discerned) are characterized by the flavor (e, ;U,r). But the neutrinos 
of definite flavor are not necessarily states of a definite mass. Instead, they are 
generally coherent superpositions of such states, 

\yd-Y.UuH)- (4) 

i 

When the standard model is extended to include neutrino mass, the mixing 
matrix U is unitary. 

In vacuum, the mass eigenstates propagate as plane waves. Leaving out the 
common phase, a beam of ultrarelativistic neutrinos with energy E at the 
distance L acquires a phase 

[^,(L))~[^,(L = 0))exp(-z^^) (5) 

Given that, the amplitude of the process —>■ is 

A(ue ^ i^e) = Uue-'^U;,, , (6) 

i 

and the probability of the flavor change for I ^ I' is the square of this ampli- 
tude. It is obvious that due to the unitarity of U there is no flavor change if all 
masses vanish or are exactly degenerate. The existence of the oscillations is a 
simple consequence of the coherence in Eq. (4). The detailed description of the 
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quantum mechanics of the oscillations in terms of wave packets is subtle and 
a bit involved (see, e.g. the reviews [22] and references therein). It includes, 
among ather things, the concept of coherence length, a distance after which 
there is no longer coherence in Eq. (4). For the purpose of this review such 
issues are irrelevant, and will not be discussed further. 

The idea of oscillations was discussed early on by Pontecorvo [23,24] and 
by Maki, Nakagawa and Sakata [25]. Hence, the mixing matrix U is often 
associated with these names and the notation Umns or Upmns is used. 

The formula for the probability is particularly simple when only two neutrino 
flavors, i^i and i^i', mix appreciably, since only one mixing angle is then relevant, 

P(i/£ — >• i^i'^i) = sin^ 26'sin^ 



1.27Am2(eV2) 



L(km) 
E,{GeY) 



(7) 



where the appropriate factors of h and c were included (the same is obtained 
when the length is in meters and energy in MeV). Here Am^ = |m| — 
(note that the sign is irrelevant) is the mass squared difference. 

Thus, the oscillations in this simple case are characterized by the oscillation 
length 

2.48E.(GeV) 
"--^^"^^ ^ Am^(eV^) ' 



and by the amplitude sin^ 29. 

For obvious reasons the oscillation studies are optimally performed at distances 
L ~ Lose from the neutrino source. At shorter distances the oscillation ampli- 
tude is reduced and at larger distances the neutrino flux is reduced making the 
experiment more difficult. Note that at distance L ^ Lose the oscillation pat- 
tern is smeared out and the oscillation probability (7) approaches (sin^26'/2) 
and becomes independent of Am^. 

The mixing matrix of 3 neutrinos is parametrized by three angles, conven- 
tionally denoted as ^12,^13,^23, one CP violating phase 5 and two Majorana 
phases ai,a2- Using c for the cosine and s for the sine, we write U as 



C12C13 S12C13 Si3e 

-■S12C23 - Ci2S23Si3e*'' C12C23 - Sl2S23Sl3e*'' S23C13 
•512523 ~ Ci2C23Si3e*'^ — C12S23 ~ •5l2C23Sl3e*'^ C23C13 



^3 



(9) 
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By convention the mixing angle du is associated with the solar neutrino os- 
cillations, hence the masses mi and m2 are separated by the smaller interval 
Am^o; (we shall assume, again by convention, that 1112 > mi) while 777.3 is 
separated from the 1,2 pair by the larger interval Am^^^^, and can be either 
lighter or heavier than mi and 7772- The situation where 7773 > 7772 is called 
'normal hierarchy', while the 'inverse hierarchy' has 7773 < 7771. Not every- 
body follows these conventions, so caution should be used when comparing 
the various results appearing in the literature. 

The general formula for the probability that the "transition" £ — > £' happens 
at L is 



i 

\m.'^-m'^\L 

=EIW/ + ^EE Wit^;,-t^€',-e^^^ ■ (10) 

i i j^i 



Clearly, the probability (10) is independent of the Majorana phases a. The 
oscillations described by the Eq.(lO) violate the individual flavor lepton num- 
bers, but conserve the total lepton number. The oscillation pattern is identical 
for Dirac or Majorana neutrinos. 

The general formula can be simplifled in several cases of practical impor- 
tance. For 3 neutrino flavors, using the empirical fact that Am'^^^^ ^ ^^Li 
and considering distances comparable to the atmospheric neutrino oscillation 
length, only three parameters are relevant in the zeroth order, the angles $23 
and and Aatm = ^^Itm^/^-^i^- However, corrections of the flrst order in 
Asoi = ^^so/-^/4-E';/ should be also considered and are included below (some 
of the terms with Agoi are further reduced by the presence of the emipirically 
smaU sin^26'i3): 

P(7/^ ^ I/^) ~ COS^ ^13 Sin^ 2^23 Sin^ Aatm (11) 

- Asol cos^ 6*13 sin^ 26123 (cos^ ^12 - sin^ 6*13 sin^ 6*12) sin 2 Aatm 

— Asol COS 5 cos ^13 sin 2^12 sin 2^13 sin 26^23 cos 26^23 sin 2 Aatm 
+ Asol sin 6 cos 6113 sin 26112 sin 26113 sin 26123 sin^ Aatm , 



Pii'i, i^e) - sin^ 2^13 sin^ ^23 sin^ Aatm ( 12) 

- Asol sin^ 6^23 sin^ 6^12 sin^ 26^13 sin 2 Aatm 

+ Asol COS 5 cos ^13 sin 2^13 sin 2^23 sin 2^12 sin 2 Aatm/ 2 

- Asol sin S cos 6I13 sin 26*12 sin 26I13 sin 26I23 sin^ Aatm , 
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P{jy^ P{yf, ^ ^e) - P{jy^ ^ j^t) , (13) 

where S is the CP phase of Eq. (9) and 
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sin^ Aatm - ^soi sin^ 6112 sin 2 A^j^ 



+ AliCos''9,ssm^29u , (14) 

where the term hnear in Ag^i is suppressed by the factor sin^ 2^13 and there- 
fore the quadratic term is also included. The disappearance probability is 
independent of the CP phase S. 

Note that the terms of first order in Agoi depend on the sign of Aatm, i-e., on 
the hierarchy. 

On the other hand, at distances much larger than the atmospheric neutrino 
oscillation length, the electron neutrino and antineutrino disappearance is gov- 
erned by 

^(i^e ^ i^x) - 1 - sin^ - cos^ ^13 [1 - sin^ 2^12 sin^ A, J , (15) 



In both of the latter cases i^x denotes any neutrino except z/g. (General expres- 
sions for mixing of three neutrinos in vacuum can be found in Eqs. (7)-(12) of 
Ref. [13].) 

In the mixing matrix (9) the presence of the phase 5 signifies the possibility 
of CP violation, the expectation that 

Pii^i' ^ i^e) + P{yi' ^ i^i) , (16) 



i.e., that for example the probability of z/^ oscillating into Vf, is different from 
the probability of z/^ oscillating into z/g- 

The magnitude of the T or CP violation is characterized by the differences 



^P{u,^Ur)-P{Ue^Pr) (17) 

= -4C13S13C23S23C12S12 sin 5[sin 2Ai2 + sin 2A23 + sin 2A31] 
= I6C13S13C23S23C12S12 sin A12 sin A23 sin A31 , 

where, as before, Aij = {ml — m|) x L/AE. Thus, the size of the effect is the 
same in all three channels, and CP violation is observable only if all three 
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masses are different (i.e., nondegenerate), and all three angles arc nonvanish- 
ing. The possibihty of CP violation in the lepton sector was first discussed in 
[26,27]. 



2.4 Neutrino oscillations in matter 

The oscillation phenomenon has its origin in the phase difference between the 
coherent components of the neutrino flavor eigenstates described by Eq. (5). 
When neutrinos propagate in matter additional contributions to the phase 
appear, besides the one caused by the nonvanishing mass of the state i/j. To 
see the origin of such phase, consider the effective Hamiltonian of neutrinos 
in presence of matter. Obviously, only phase differences are of importance. 

Without invoking any nonstandard interactions two effects are present. All 
active neutrinos interact with quarks and electrons by the neutral current 
weak interactions {Z exchange), but only electron neutrinos and antineutrinos 
interact with electrons by the exchange of W. The corresponding effective 
potential is ^ {Vc stands for the charged current) 

Vc{iye) = V2GFN,, and Vc(z>e) = -x/SG'i.iVe , (18) 

where Ng is the electron number density. (There are no or r leptons in normal 
matter, hence there is no analogous potential for i/f^ or 1/-,-.) In practical units 

Similarly, any active neutrino acquires an effective potential due to the neutral- 
current interaction Vn = —GpN^/y/^. The corresponding effective potential 
is absent for sterile neutrinos. 

Given the effective potential, Eq.(18), electron neutrinos travelling distance L 
in matter of constant density A^e acquire an additional phase 

i^eiL) = t/e(0)e-^^^^^^^ . (20) 



^ The effective potential was introduced first by Wolfcnstcin [28], and used also in 
Refs. [29,30] who corrected the missing ^/2 in the original paper. Finally, the correct 
sign was obtained in Ref.[31]. 
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(19) 



The corresponding matter oscillation length [28] is therefore 

^ _ 27r _ 1-7 X lO^(m) 
° ~ ^/2GfN, ~ p(g/cm3) n ■ 



(21) 



Unlike the vacuum oscillation length, Eq.(8), the matter oscillation length Lq 
is independent of the neutrino energy. Note that the matter oscillation length 
in rock is Lq ~ 10^ km, and in the center of the Sun Lq ~ 200 km. 

Considering for simplicity just two mass eigenstates vi and V2 that are com- 
ponents of the flavor eigenstates and z/q, with the mixing angle 6, we obtain 
the time (or space) development Schrodinger equation 




Vcsc 



\ 



Vcsc 



m| 
2E 



^ + VcS^ 



1 









(22) 



where, as before c = cos^ and s — sin. 6. 

The 2x2 matrix above can be brought to the diagonal form by the transfor- 
mation 



^Im = COS 0m " Siu 9^. 
T^2m = T^e Siu 9m + COS 9m , 



(23) 



where the new mixing angle in matter, 9m, depends on the vacuum mixing 
angle 9 and on the vacuum and matter oscillation lengths Lose and Lq, 



tan 2^™ =tan2^ 1 



Lq cos 29 



(24) 



The effective oscillation length in matter is then 



sin 29m 
' sin 2^ 



1 + 



2L, 



cos29 



-1/2 



(25) 



and the probability of detecting z/g at a distance L from the v^. source has the 
usual form, but with 9 ^ 9m and L^sc Lm, 



P{E^, L, 9, Am^) = 1 - sin^ 29m sin^ — . 



(26) 



When considering oscillations with two flavors, the mixing angle 9 can be re- 
stricted to the interval (0, tt/2). In vacuum, only sin^ 29 is relevant, and hence 



12 



only half of that interval, (0, 7r/4) could be used. However, once matter oscilla- 
tions are present, also the cos 29 becomes relevant, and thus the whole (0, 7r/2) 
interval might be needed. The part of the parameter space corresponding to 
the (7r/4,7r/2) angles was called the "dark" side in Ref.[32], who suggested 
using tan^ 9 in the plots instead of sin^ 29. (The advantage of tan^ 9 is that 
when plotted on the log scale, the reflection symmetry around tan^ 9 = 1 for 
vacuum oscillations is maintained.) The importance of allowing the full range 
of the mixing angles in the three flavor analysis was stressed earlier, see e.g. 
Ref. [33]. 

With our convention m2 > rrii mixing angles in the interval (7r/4,7r/2) would 
mean that contains dominantly the heavier component V2- (This is not the 
case in practice, see Sect. 3.1.) It appears now that at least two of the three 
mixing angle are large (and < 7r/4), and the constraint on the third one, ^13 
are still such that hnear plots are more reveahng. Therefore, in the following 
we use the traditional plots with sin^ 29. 

The same results can be obtained, naturally, by rewriting the equation of 
motion in the flavor basis 

, / , . \ / 008 261 I 1 sin 26* \ / , , \ 

Here one can see clearly that in matter, unlike in vacuum, the oscillation 
pattern depends on whether the mixing angle 9 is smaller or larger than 7r/4. 
(For antineutrinos the sign in front of 1/Lq is reversed.) 

In matter of a constant density we can now consider several special cases: 

• Low density limit. Lose Lq. In this case matter has a rather small effect. 
On Earth, one is able to observe oscillations only provided Lose < Earth di- 
ameter, and therefore this limit applies. The matter effects could be perhaps 
observed as small day-night variation in the solar neutrino signal. In long 
baseline oscillation experiments matter effects could cause difference in the 
oscillation probability of neutrinos and antineutrinos, hence these effects, 
together with the hierarchy problem (i.e., whether > m2 or 777,3 < it^i) 
arc crucial in the search for CP violation and its interpretation. 

• High density limit. Lose ^ Lq. The oscillation amplitude is suppressed by 
the factor Lo/\Losc\- For 777^ > ml the matter mixing angle is 9jn — > 90° and 
thus — > ^2- 

• I Lose I ~ -^0- In this case the matter effects can be enhanced. In particular, 
for Lose/Lo = cos2^^ one has sin^ 2^^^ = 1, i.e. the maximum mixing, even 
for small vacuum mixing angle 9. This is the basis of the Mikheyev-Smirnov- 
Wolfenstein (MSW) effect [28,34]. 
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When neutrinos propagate in matter of varying density, the equations of mo- 
tion (22) or (27) must be solved. There is a vast literature on the subject, in 
particular on the application to the solar neutrinos. There, the v^. are often 
created at densities above the resonance density, and therefore are dominantly 
in the higher mass eigenstate. When the resonance density is reached, the two 
instantaneous eigenvalues are almost degenerate. In the adiabatic regime, the 
neutrino will remain in the upper eigenstate, and when it reaches vacuum it 
will be (for small vacuum mixing angle &) dominantly in the state that we 
denoted above as v^- In the general case, there is a finite probability for 
jumping from one eigenstate to the other one, and the conversion might be 
incomplete. The average survival probability is [35] 

(P(l/e ^ Ve)) = ^[1 + (1 - 2P^) COs2^„(p„„^) COS 2^] , (28) 



Usually C0s2^TO(Pmaa:) — "1 and thus {P{ye — i^e)) — Sin^ Q + Pa; COS 2^. 

The transition point between the regime of vacuum and matter oscillations is 
determined by the ratio 



2\/2GFNeE^ 



0.22 



1 MeV 



lOOg/cm^ 



7 X lO-^eV^ 



(29) 



If this fraction is larger than unity, the matter oscillations dominate, and when 
this ratio is less than cos 29 the vacuum oscillations dominate. Generally, there 
is a smooth transition in between these two regimes. 

The electron neutrino survival probability is illustrated in Fig. 1, where it 
is plotted against iE^/Am^, see Eq.(29). In Fig. 1 it was assumed that the 
neutrinos originated in the center of the Sun, hence the relatively sharp feature 
at 4:E^/Am'^ ~ 2.6 x 10^^ eV"^ where according to the Eq.((29) Losc/Lo = 1 
for the central density of the Sun. Note that below and above this dividing 
line the survival probability is almost independent of the neutrino energy, 
hence no spectrum distortion is expected. The subsequent increase in the 
neutrino survival probabihty for larger values of AE^/Am? is caused by the 
'nonadiabatic' transition, i.e. by the gradual increase of the jump probability 
towards the limiting value Pj; — > 1. 

For the parameters corresponding to the preferred solar solution (sin^ 29 ~ 
0.83 and Am^ ~ 7 x 10"^ eV^) the pp neutrinos with AE^/Am"^ < 2.4 x 10^° 
eV^^ and the '^Be neutrinos with AEj,/ Am? = 4.9 x lO-*^*^ eV~^ undergo vacuum 
oscillations, while the neutrinos with AEy/ Arn^ > 2.55 x 10^^ eV~^ undergo 
MSW matter oscillations. (Clearly, that must be the case since otherwise it 
would be impossible to understand the ~ 0.3 suppression of the ^B neutrino 
i/g fiux observed in the charged current reactions, see 3.1.) 
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4EyAm^ (eV"') 



Fig. 1. Schematic illustration of the survival probability of v^. created at the solar 
center. The curves are labelled by the sin^ 29 values. 



2. 5 Tests of CP, T and CPT invariance 



In vacuum CP conservation implies that P(z/^ Ufj) = P{ue ue') (see 
Eq.(18). Violation of this inequality thus would mean that CP is not conserved 
in the lepton sector. Substantial effort is devoted to the CP tests. Matter 
effects can induce inequality between P^ui z/^') and P{u£ z/^) and so the 
analysis must carefully account for them. 

If CP is not conserved, but CPT invariance holds, then T invariance will be 
also violated, i.e. P(z/£ ugi) — P{v(^i — >■ z/^) need not vanish. Here matter 
effects cannot mimic the apparent T invariance violation. 

Finally, if CPT is not conserved, then P{v(, Ufi) — P{jye' —* ^i) might be 
nonvanishing. Many tests of the CPT invariance in the neutrino sector have 
been suggested, see e.g. [36] or references listed in [13]. 

CPT invariance is based on Lorentz invariance, hermiticity of the Hamiltonian 
and locality. Its violation would have, naturally, enormous consequences. Yet 
there are many proposed scenarios of CPT violation, in particular in the 
neutrino sector (for a whole series of papers on that topic see e.g. [37]). 

CPT invariance implies that neutrino and antineutrino masses are equal. If 
that is not true, then the Am? as determined in the solar neutrino experiments 
(thus involving z/g) might not be the same as the Am^ needed to explain the 
LSND result which involves z/^ z/g oscillations, see 3.1. 4^ below. That was the 
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gist of the phenomenological proposal in Ref . [38] . (See also further elaboration 
in [39].) 



With the demonstration of the consistency between the observed solar v^. 
deficit and the disappearance of reactor Pg by the KamLAND collaboration 
(see 5. i. 5 below), this possibility seems unlikely, even though a proposal has 
been made to accommodate CPT violation in that context [40], see also [42]. 
As has been shown in [41], the consistency of the solar v oscillation solution 
and the KamLAND reactor result can be interpreted as a test of CPT giving 



where Am^ and Am? refer to the mass cigenstates i^i and 1^2 involved in the 
observed solar and reactor neutrino oscillations. 

To test for the CP invariancc violation experimentally, one would compare the 
probabilities P(z/^ Ve) and ^(2?^ Ve)- This could be done realistically with 
~ 1 GeV Ufj_ beams at a distance L ~ Ei./Am'^^^ such that the contribution 
involving Am^^,; are small. The effect of matter, however, must be included. 
Using the notation 

Sij = sin%, Cij = cosOij, Aij = AmfjL/AEi, (31) 
we obtain the formula 



Here the first term gives the largest effect, while the terms in the third and 
fourth line represent small CP conserving corrections (proportional to sin A21 
and sin^ A21. The term with sin 5 in the second line violates CP symmetry, 
while the term with cos 5 preserves it. Finally, the term with aL/AEn in the 
last line represents the matter effects. 

The matter effects are characterized by 



a = 2V2GfN^E^ = 1.54 x 10-^FeP(g/cm=^)^^(GeV) (a is in (eV^)). (33) 



Ami - ^^ll < 1-3 X 10-2 eV^ 90% CL , 



(30) 




•31 J 



(32) 



The probability -P(t'^ 



i/g) is obtained by the substitution 6 



—6 and 



a 



—a. 
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To test for CP symmetry, one would determine 



Acp 




sin 2^12 . J. aL cos A; 



■32 



(34) 



•31 



where we used the empirical fact that cos ^23 ^ sin ^23 and sin(A2i) ~ A21 
for the distances and energies usually considered. Since ^^13 is small, the CP 
asymmetry can be enhanced. However, the individual terms in Eq. (34) depend 
on 6*13, so for smaller 6*13 it is more difficult to reach the required statistical 
precision. 

There are several parameter degeneracies in Eq. (34) when separate measure- 
ments of P(t'^ —>■ I'e) and Pii'lj, Vf.) are made at given L and E^. {i) There 
can be two values 5, ^13 and 5', 6'^^ leading to the same probabilities, {ii) Sign 
of A31 and A32 (i.e. the normal or inverted hierarchy), where one set 5, ^13 
gives the same oscillation probabilities with one sign, as another set 5' , O'^.^ 
with the opposite sign of A^i. (in) Since the mixing angle ^23 is determined 
in experiments sensitive only to sin^ 2^23 there is an ambiguity between ^23 
and 7r/2 — ^23- However, for the preferred value ^23 ~ 7r/4 this ambiguity is 
essentially irrelevant. Various strategies to overcome the parameter degenera- 
cies have been proposed. The choice of the neutrino energy E,, (and whether a 
wide or narrow beam is used) and the distance L play essential role. Clearly, 
if some of the so far unknown parameters (e.g. ^13 or the sign of the hier- 
archy) could be determined independently, some of these ambiguities would 
be diminished. (For some of the suggestions how to overcome the parameter 
degeneracy see e.g. [43,44].) 

2. 6 Violation of the total lepton number conservation 

Neutrino oscillations described so far are insensitive to the transformation 
properties under charge conjugation, i.e., whether neutrinos are Dirac or Ma- 
jorana particles. However, if the mass eigenstates Vi are Majorana particles, 
then z/ — s> z/ oscillations that violate the total lepton number conservation are 
possible. There are two kinds of such processes. 

Recall that in the standard Model charged current processes the neutrino 
produces the negatively charged lepton £" while antineutrino produces 
the positively charged antilepton The standard model also requires that 
in order to produce the lepton the neutrino z/£ must have (almost) purely 
negative helicity, and similarly for the P^ must have (almost) purely positive 
helicity. The amplitudes for the 'wrong', i.e., suppressed helicity component 
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is only of the order niu/Ey and therefore vanishes for masslcss neutrinos. But 
we know now that neutrinos are massive, although light, and the therefore 
'wrong' helicity states are present. 

When neutrinos are Majorana particles, the total lepton number may not be 
conserved. Thus, neutrinos vg, born together with the leptons can create 
Icptons (or even a different flavor as long as the helicity rules are 
obeyed. The amplitude of this process is of the order rriy/Ey, thus small, 
independent of the distance the neutrinos travel. Therefore, this first kind of 
the 1/ — > P transformation should not be really called oscillations. When such a 
transformation happens inside a nucleus, it leads to the process of neutrinoless 
double beta decay discussed in more detail below. 

For Majorana neutrinos, there are several nonstandard processes involving he- 
licity fiip in which left-handed neutrinos vj^ are converted into right-handed 
(anti) neutrinos z/^. This can happen for Majorana neutrinos with a transition 
magnetic moment /ij-,-. In a transverse magnetic field B±, can be connected 
to u'^ji- However, the transition probability is proportional to the small quan- 
tity \iivB\ that vanishes for massless neutrinos. There are also models that 
predict neutrino decay involving a helicity-fiip and a Majoron x production. 
Again, the decay rate is expected to be small. If such processes exist one ex- 
pects, among other things, that solar z/g could be subdominantly converted 
into Ve- The recent limit on the solar z/g fiux, expressed as a fraction of the 
Standard Solar Model ®B fiux is 2.8 x 10"^ [45]. (See also [45] for references 
to some of the theoretical model expectations for such — > P conversion.) 

In addition, there could be also transitions without helicity fiip, which require 
that both Dirac and Majorana mass terms arc present. These second class os- 
cillations [46] involve transitions z/^ — > z/£, i.e., the final neutrino is sterile and 
therefore unobservable. Estimates of the experimental observation possibilities 
of the neutrino-antineutrino oscillations are not encouraging [47] . 

Study of the neutrinoless double beta decay (0z//9/9) appears to be the best 
way to establish the Majorana nature of the neutrino, and at the same time 
gain valuable information about the absolute scale of the neutrino masses. 
Double beta decay is a rare transition between two nuclei with the same mass 
number A involving change of the nuclear charge Z by two units. The decay 
can proceed only if the initial nucleus is less bound than the final one, and 
both must be more bound than the intermediate nucleus. These conditions are 
fulfilled in nature for many even-even nuclei, and only for them. Typically, the 
decay can proceed from the ground state (spin and parity always 0"*") of the 
initial nucleus to the ground state (also 0"'') of the final nucleus, although the 
decay into excited states (0+ or 2+) is in some cases also energetically possible. 
In Fig. 2 we show a typical situation for A = 136. There are 11 candidate nuclei 
(all for the /3~P~ decay) with Q value above 2 MeV, thus potentially useful 
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for the study of the 0uj3f3 decay. (Large Q values are preferable since the rate 
scales like and the background suppression is typically easier for larger Q.) 




A=136 




136 
59 



Pr 



- 7 



- 6 



- 5 
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■57 



La 13.6, 

58 



Ce 



(MeV) 



136 
56 



Ba 



Fig. 2. Masses of nuclei with A = 136. The even-even and odd-odd nuclei are con- 
nected by dotted lines. ^^^Xe is stable against ordinary /3 decay, but unstable against 
(3~ (3~ decay. The same is true for ^^^Ce, however, the decay is expected to 

be slower than the /3^/3^ decay. 

There are two modes of the double beta decay. The two- neutrino decay, 2i//3/3, 



conserves not only electric charge but also lepton number. On the other hand, 
the neutrinoless decay. 



violates lepton number conservation. One can distinguish the two decay modes 
by the shape of the electron sum energy spectra, which are determined by the 
phase space of the outgoing hght particles. Since the nuclear masses are so 
much larger than the decay Q value, the nuclear recoil energy is negligible, 
and the electron sum energy of the Ou^P is simply a peak at T^i + Tg2 — Q 
smeared only by the detector resolution. 

The 2i//?/? decay is an allowed process with a very long lifetime ~ 10^° years. 
It has been observed now in a number of cases [12]. Observing the 2i'f3f3 decay 
is important not only as a proof that the necessary background suppression 
has been achieved, but also allows one to constrain the nuclear models needed 
to evaluate the corresponding nuclear matrix elements. 



{Z, A)^{Z + 2, A) + + + Pel + i>e2 



(35) 



{Z,A) (Z + 2,A) + er + 62 



(36) 
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The Oz//3/3 decay involves a vertex changing two neutrons into two protons 
with the emission of two electrons and nothing else. One can visualize it by 
assuming that the process involves the exchange of various virtual particles, 
e.g. light or heavy Majorana neutrinos, right-handed current mediating Wr 
boson, SUSY particles, etc. No matter what the vertex is, the OuPP decay 
can proceed only when neutrinos are massive Majorana particles [48]. In the 
following we concentrate on the case when the Ou^jS decay is mediated by 
the exchange of light Majorana neutrinos interacting through the left-handed 
V — A weak currents. The decay rate is then. 



0+)]-i = G°^(Eo,^) 



9a 



(37) 



where G'^'^ is the exactly calculable phase space integral, {mpp) is the effective 
neutrino mass and MqIj,^ Mp" are the nuclear matrix elements. 

The effective neutrino mass is 

{mpp)^\T.\Uei\'m,,e'^^\ , (38) 



where the sum is only over light neutrinos (mj < 10 MeV) . The Majorana 
phases ccj were defined earlier in Eq.(9). If the neutrinos Ui are CP eigenstates, 
ai is either or tt. Due to the presence of these unknown phases, cancellation 
of terms in the sum in Eq.(38) is possible, and {mpp) could be smaller than 
any of the ml,^ . 

The nuclear matrix elements, Gamow-Teller and Fermi, appear in the combi- 
nation 

M^T - = (/I T.H{rm EmWr^ U ■ - 4k)) • (39) 

9a Ik \ 9a J 



The summation is over all nucleons, |i),(|/)) are the initial (final) nuclear 
states, and H(rik, Em) is the 'neutrino potential' (Fourier transform of the 
neutrino propagator) that depends (essentially as 1/r) on the internucleon 
distance. When evaluating these matrix elements the short-range nucleon- 
nucleon repulsion must be taken into account due to the mild emphasis on 
small nucleon separations. 

There is a vast literature devoted to the evaluation of these nuclear matrix 
elements, going back several decades. It is beyond the scope of the present re- 

^ The same quantity is sometimes denoted as (m,^) or niee 



20 



view to describe this effort in detail. The interested reader can consult various 
reviews on the subject, e.g. [12,49,50]. 

Obviously, any uncertainty in the nuclear matrix elements is reflected as a 
corresponding uncertainty in the (m^^). There is. at present, no model inde- 
pendent way to estimate the uncertainty, and to check; which of the calculated 
values are correct. Good agreement with the k;nown 2i'j3f3 is a necessary hut 
insufficient condition. The usual guess of the uncertainty is the spread, by a 
factor of ~ 3, of the matrix elements calculated by different authors. Clearly, 
more reliable evaluation of the nuclear matrix element is a matter of consider- 
able importance. (For a recent attempt to reduce and understand the spread 
of the calculated values, see [51].) 

The OuPP decay is not the only possible observable manifestation of lepton 
number violation. Muon-positron conversion, 

fi- + {A, Z)^e+ + {A, Z-2) , (40) 

or rare kaon decays K^^^,, Keen and K^^tt, 

K+ ^+^+7r- , K+ e+e+TT" , K+ /i+e+Tr" , (41) 

are examples of processes that violate total lepton number conservation and 
where good limits on the corresponding branching ratios exist. (See Ref. [52] 
for a more complete discussion.) However, it appears that the Qvf3f3 decay is, 
at present, the most sensitive tool for the study of the Majorana nature of 
neutrinos. 

2. 7 Direct measurement of neutrino mass 

Conceptually the simplest way to explore the neutrino mass is to determine 
its effect on the momenta and energies of charged particles emitted in weak 
decays. 

In two-body decays, e.g. tt"*" — > ^'^ + the analysis is particularly straightfor- 
ward, at least in principle. In the system where the decaying pion is at rest, 
the energy and momentum conservation requirements mean that 

ml = ml + ml- 2m^^ml+pl ■ (^2) 

However, the neutrino mass squared appears as a difference of two very large 
numbers. Hence the uncertainties in tUt^^p^ and even mean that the cor- 
responding mass limit is only rUi, < 170 keV. 
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This problem can be avoided by studying the three body decays, in particular 
the nuclear beta decay. Near the endpoint of the beta spectrum a massive 
neutrino has so little kinetic energy that the effects of its finite mass become 
more visible. The electron spectrum of an allowed beta decay is given by the 
corresponding phase-space factor 

^ ~ F{Z, E,)p,E,{Eo - EMEo - E^f - ml]'/' , (43) 

where E^iPe is the electron energy and momentum and F{Z,Ee) describes 
the Coulomb effect on the outgoing electron. The quantity Eq is the endpoint 
energy, the difference of total energies of the initial and final systems. Clearly, 
the effect of finite neutrino mass becomes visible if {Eq — Eg) ~ m^, i.e. very 
near the threshold. 

For the case of several massive neutrinos with mixing, the beta decay spectrum 
is an incoherent superposition of spectra like (43) with corresponding weights 
|^7eiP for each mass eigenstate ml.. If the experiment has insufficient energy 
resolution the quantity 

mlf^-l:\Ue^\'ml^ (44) 

i 

(using the RPP notation [20]) could be determined from the electron spectrum 
near its endpoint, where the sum is over all the experimentally unresolved 
neutrino masses m^,.. 

Based on an upper limit for the ml^^^^ one can deduce several limits that do 
not depend on the mixing parameters |C/ei|^- First, at least one of the neutrinos 
(i.e. the one with the smallest mass) has a mass less or equal to that limit, 
ml^.^ < ml^^^\ Moreover, if all (with an emphasis on all) \Aml-\ values are 
known, an upper limit of all neutrino masses is mf, < mfj-'^^^ + y],v-o I Am? |. 
Thus, if we assume that the |Am^j | values deduced from the experiments 
on solar (and reactor) and atmospheric oscillation studies (and include also 
the LSND result) cover all possibihties, and combine that knowledge with the 
j^2(eff) from tritium beta decay, we may conclude that no active neutrinos 
with mass more than 'a few' eV exists. 



3 Experimental Results and Interpretation 

Positive evidence for neutrino mass has so far been obtained only in mea- 
surements of neutrino oscillations. As noted above in 2.3, the neutrino masses 
enter only through the differences in squared masses (i.e.. Am?- = ml — mf), 
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so although these measurements provide lower limits to the mass values (e.g., 
mf > Amjj) they do not actually determine the masses. As discussed in 2.6 
and 2.7, direct neutrino mass measurements and neutrinoless double beta 
decay allow determination of quantities related to the values of the masses 
themselves (as opposed to differences). However, these experiments have so 
far yielded only limits. Finally, additional constraints have been derived from 
measurements of the cosmic microwave background radiation and galaxy dis- 
tribution surveys through the effect of neutrinos on the distribution of matter 
in the early universe. In this section, we summarize the positive observations 
from oscillation measurements followed by brief discussion of the direct mea- 
surements, double beta decay searches, and recent constraints from cosmology. 



3.1 Neutrino Oscillation Results 

The first hints that neutrino oscillations actually occur were serendipitously 
obtained through early studies of solar neutrinos and neutrinos produced in 
the atmosphere by cosmic rays ("atmospheric neutrinos"). In fact, the at- 
mospheric neutrino measurements were a byproduct of the search for proton 
decay using large water Cerenkov detectors. So it is somewhat ironic that 
although there was substantial interest in searching for neutrino oscillations, 
the first evidence for this phenomena came from experiments designed for very 
different purposes. 

As shown below, recent studies definitively establish that the solar neutrino 
flux is reduced due to flavor oscillations, and so it is now clear that the first real 
signal of neutrino oscillations was the long-standing deficit of solar neutrinos 
observed by Ray Davis and collaborators using the Chlorine radiochemical 
experiment in the Homestake mine. While it took almost three decades to 
demonstrate the real origin of this deficit, the persistent observations by Davis 
et al. and many other subsequent solar-z/ experiments were actually indications 
of neutrino oscillations. We will discuss this subject in more detail below. 

3.1.1 Atmospheric Neutrinos 

The Kamiokande experiment in Japan [53] and the 1MB experiment in the 
US [54]were pioneering experimental projects to develop large volume water 
Cerenkov detectors with the primary goal of detecting nucleon decay, as pre- 
dicted by Grand Unified Theories developed in the 1970's [55]. Although these 
detectors were located deep underground to avoid cosmic ray-induced back- 
ground, they both encountered the potential background events produced by 
atmospheric neutrinos (both and u^) that easily penetrated to these sub- 
terranean labs and (rarely) produced energetic events in the huge detectors. 
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Fig. 3. Distribution of observed atmospheric neutrino events vs. zenith angle from 
the SuperKamiokande experiment, compared with Monte Carlo simulations [60]. 
The blue hatched region represents the prediction without neutrino oscillations and 
the red line includes the effect neutrino oscillations. (PC means 'partially contained' 
events.) 



And indeed, both experiments [56,57,58] (along with the Soudan experiment 
[59]) observed that the ratio of i/^-induced events to z/g-induced events was 
substantially reduced from the expected value of ~ 2. The decay chain of tt^ 
produced in the upper atmosphere would produce (through the subsequent 
/i-decay) a u^, z/^, and a z/g (or z/g). Thus, based on rather simple basic ar- 
guments one expects the ratio of z/^/z/g events to be about ~ 2 - and this 
is supported by more detailed Monte Carlo simulations. The observed values 
were closer to ~ 1, which was viewed as an anomaly for many years. Here 
again, although z/-oscillations could clearly cause this anomaly there was not 
enough corroborative evidence to substantiate this explanation. 



However, the situation dramatically changed in 1998 when the larger exper- 
iment, Super-Kamiokande, reported a clearly anomalous zenith angle depen- 
dence of the z/p events [61]. The measurements indicated a deficit of upward- 
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going z/^-induccd events (produced ~ 10^ km away at the opposite side of 
the earth) relative to the downward-going events (produced ~ 20 km above). 
The Vf. events displayed a normal zenith angle behavior consistent with Monte 
Carlo simulations. Since at that time (1998) the solar- i/ problem was still un- 
resolved, these data represent the first really sohd evidence for i/-oscillations. 
More recent measurements [60] are displayed in Fig. 3, and the conclusion that 
flavor oscillations are responsible is essentially inescapable. Moreover, the de- 
duced values of sin^ 26* > 0.90 (90% C.L.) indicate a surprisingly strong mixing 
scenario where the muon-type neutrino seems to be a fully- mixed superposi- 
tion of all three mass eigenstates. (This situation is completely contrary to 
the quark sector, where the mixing between generations is generally small.) 
The most recently reported value of Am^ derived from the Super-Kamiokande 
results is 0.0020 eV^ [62] (for the error bars, see Table 1). 

The observation of the angular distribution of upward-going muons produced 
by atmospheric neutrinos in the rock below the MACRO detector [63] supports 
the conclusion that the observed effect is due to the — > oscillations, and 
disfavors — > v sterile assignment. 

Although the parallel effort to detect neutrino disappearance at nuclear reac- 
tors had made steady progress (setting upper limits and establishing exclusion 
plots) for many years, these experiments made a strong contribution at this 
point: the failure to observe Vf. disappearance at CHOOZ [64] and Palo Verde 
[65] in the region near Am^ ~ 0.0020 eV^ imphes that the i/^ disappearance 
observed by Super-Kamiokande does not involve substantial appearance. 
(This inference assumes that my = rrip for each eigenstate as required by 
CPT invariance.) Thus, it would seem that the z/^'s must be oscillating into 
or other more exotic particles, such as "sterile" neutrinos. 



3.1.2 Solar Neutrinos 

The interpretation of solar neutrino measurements involves substantial input 
from solar physics and the nuclear physics involved in the complex chain of 
reactions that together are termed the "Standard Solar Model" (SSM) [66] . 
The predicted flux of solar neutrinos from the SSM is shown in Fig. 4 as a 
function of neutrino energy. The low energy p-p neutrinos are the most abun- 
dant and, since they arise from reactions that are responsible for most of the 
energy output of the sun, the predicted flux of these neutrinos is constrained 
very precisely (±2%) by the solar luminosity. The higher energy neutrinos are 
more accessible experimentally, but the fluxes are less certain due to uncer- 
tainties in the nuclear and solar physics required to compute them. 

The early measurements included radiochemical experiments sensitive to in- 
tegrated Ug flux such as the Chlorine [68] (threshold energy 0.814 Mev) and 
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Neutrino Energy (MeV) 



Fig. 4. Energy distribution of the flux of solar neutrinos predicted by the standard 
solar model [66], as computed by J. Bahcall [67]. The ranges of energies associated 
with the various experiments are indicated at the top of the figure. 

Gallium [69] (0.233 MeV) experiments. Live counting was developed by the 
Kamiokande [70] and then the SuperKamiokande [71] experiments, based on 
neutrino-electron scattering, enabling measurements of both the flux and en- 
ergy spectrum. Together these experiments sampled the solar neutrino flux 
over a wide range of energies. As can be seen in Fig. 8 below, all these ex- 
periments reported a substantial deficit in neutrino flux relative to the SSM. 
While it was realized that neutrino oscillations could be an attractive solution 
to this problem, it was problematic to establish this explanation with certainty 
due to the dependence on the SSM, its assumptions, and sensitivity to input 
from nuclear and solar physics. 

However, it is now clear that the precise measurements of solar neutrino fluxes 
over a range of energies coupled with the amazing flavor transformation prop- 
erties of neutrinos validates the SSM in a beautiful and satisfying manner. 
The solar-i/ experiments and the SSM have been reviewed in detail elsewhere 
[20], and so it is not appropriate to repeat the details here. We present these 
solar neutrino flux measurements together in Fig. 8, where we also display 
the excellent agreement with the best fit solution to the combined solar- z/ and 
KamLAND reactor data [72] to demonstrate the remarkably coherent picture 
that is strongly supported by these impressive experimental measurements 
and the associated SSM input. 
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3.1.3 SNO and Kami AND 



With the advent of the new millenium, the stage was set for a synthesis of 
the study of solar neutrinos using a powerful new water Cerenkov detector 
(Sudbury Neutrino Observatory, SNO) with the study of the disappearance 
of reactor antineutrinos using a large scintillation detector located deep un- 
derground (KamLAND). The results of these experiments provide definitive 
evidence that the solar deficit is indeed due to fiavor oscillations, and that this 
effect is demonstrable in a "laboratory" experiment on earth. 

The SNO experiment combines the now high-developed capability of water 
Cerenkov detectors with the unique opportunities afforded by using deuterium 
to detect the solar neutrinos [73,74]. Low energy neutrinos can dissociate deu- 
terium via the charged current (CC) reaction 

v^ + d^ e~ +2p (45) 

or the neutral current (NC) reaction 

ve,-\-d^Vi-\-p-\-n. (46) 



Only i/e can produce the CC reaction, but all flavors £ — e, Ijl,t can contribute 
to the NC rate. The CC reaction is detected via the energetic spectrum of e~ 
which closely follows the solar z/g spectrum. The NC reaction involves three 
methods for detection of the produced neutron: (a) capture on deuterium and 
detection of the 6.25 MeV 7-ray, (b) capture on CI (due to salt added to the 
D2O) and detection of the 8.6 MeV 7-ray, or (c) capture in ^He proportional 
counters immersed in the detector. There are also some events associated with 
the elastic scattering of the solar-z/ on e~ in the detector which is dominated by 
the charged current reaction (again only z/g) but has some ~ 20% contribution 
from neutral currents (all flavors equally contribute). 

The SNO collaboration has published data on the CC and NC rate (from pro- 
cesses (a) and (b)). Additional data from the NC process (c) will be forthcom- 
ing in the future. Nevertheless, the reported results (see Fig. 5) demonstrate 
very clearly that the total neutrino flux (i/g + i^// + i^r as determined from NC) 
is in good agreement with the SSM. but that the z/g flux is suppressed (as 
determined from CC). This represents rather definitive evidence that the z/g 
suppression is due to fiavor-changing processes that convert the z/g to the other 
flavors, as expected from z/-oscillations. Furthermore, the observed value of i/g 
flux and the observed energy spectrum, when combined with the other solar-i/ 
measurements strongly favor another large mixing angle scenario at a lower 
value of Am^ ~ 10~^ eV^. 

The KamLAND experiment [75] represents a major advance in the develop- 
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Fig. 5. Measured solar neutrino fluxes from SNO for the NC and CC processes [74], 
along with elastic scattering events (ES) and the SSM prediction [66]. 

ment of reactor z/g measurements. In order to reach the low values of Am^ ~ 
10~^ eV^ indicated by the solar-z/ data, distances of order ~ 200 km are nec- 
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KamLAND experiment [75]. The shaded region is that expected due to neutrino 
oscillations with large mixing angle parameters as determined from solar neutrino 
data. 
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Fig. 7. Region of parameter space constrained by simultaneous fit to solar-z^ and 
KamLAND data, from [72]. The best fit values are Am^ = 7.1 X 10-5 eV^ and 
tan2(9 = 0.41. 

essary (given the fixed range of Pe energies from reactors). The loss of rate 
due to scaling is severe, which requires substantial increases in source 
strength (i.e., reactor power) and detector size. Such a huge detector, sensi- 
tive to low-energy inverse beta-decay events from reactor P^, would be very 
susceptible to background from cosmic radiation and so must be located deep 
underground. By fortunate coincidence, the old Kamiokande site in Japan is 
very deep (~ 1000 mwe) and located an average distance of ~ 200 km from 
a substantial number of large nuclear power reactors. Thus the KamLAND 
experiment, a large liquid scintillator detector, was built at this site to study 
the disappearance of Ue from nuclear reactors. For the first time in the long 
history of reactor Ue experiments (dating back to the original discovery of the 
Ue by Reines and Cowan [2]) a substantial deficit in event rate was observed 
(Fig. 6). In fact this deficit is just as predicted by the solar solution to the 
solar-z/ oscillation solution, and the more precisely constrained values of Am^ 
and sin^ 29 from a global analysis of the solar-i/ and KamLAND data are 
shown in Fig. 7. 

As mentioned previously. Fig. 8 shows a summary of the solar-i/ data compared 
with the SSM with and without neutrino oscillations. The plotted experiments 
are: Ga, combined Gallium measurements from GALLEX and SAGE [69]; CI, 
Chlorine measurement from Homestake mine [68]; SNO, Sudbury Neutrino 
Observatory (CC and NC) [72,73,74]; and SK, Super-Kamiokande[71]. In this 
plot one can clearly see the decreasing survival fraction with increasing energy 
in the progression Ga— >^Cl^SNOcc (all sensitive only to the component). 
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The SNOnc measurement shows no suppression, whereas the SK data exhibit 
the intermediate suppression due to the partial contribution of NC events to 
their elastic scattering signal. 

In summary, the experimental studies of solar neutrinos, atmospheric neutri- 
nos, and reactor antineutrinos have established neutrino oscillations with two 
different mass scales, Am^ ~ 2.0 x 10~^ eV^ and Am^ ~ 7.1 x 10~^ eV^, both 
with large associated mixing angles. The allowed regions are shown in Fig. 9, 
and together these experiments constrain many of the matrix elements in the 
3x3 mixing matrix for the neutrinos along with the mass differences Amfg 
and Am32. The results for these parameters are listed in Table 1. 




Fig. 8. Ratio of solar neutrino flux to SSM (without neutrino oscillations) for vari- 
ous experiments (see text). Filled circles are experimental data (with experimental 
uncertainties only) and open circles are theoretical expectations based on SSM with 
best fit parameters to KamLAND and solar-z/ data (uncertainties from SSM and 
oscillation fit combined). All charged current experiments show a substantial deficit 
and all are in excellent agreement with the expected values. 



Table 1 

Neutrino Oscillation Parameters Determined From Various Experiments (2003) 



Parameter 



Value iblo" 



Reference Comment 



C12 

sin2 2^23 
sin^ 26*13 



eV2 



7-lia6 X 10- 

on co+2.4 
"^^•"^-2.3 



2.0toA X 10"^ eV2 



> 0.94 
< 0.11 
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[72] 
[62] 
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Fig. 9. Allowed regions of parameter space (90% C.L.) determined by atmospheric 
neutrino measurements [62] (2^3 mixing) and by solar-z/ and reactor-z^ measure- 
ments [72] (1 <-> 2 mixing). 

3.1.4 LSND 

There is one other experiment that claims to observe neutrino oscillations: 
the Liquid Scintillator Neutrino Detector (LSND) [76] at Los Alamos. In this 
experiment, the neutrino source was the beam dump of an intense 800 MeV 
proton beam where a large number of charged pions were created and stopped. 
Since the vr" capture on nuclei with very high probability, essentially only the 
71+ decay, producing and then and z/g (from /i"*" decay). These neutrinos 
all have very well defined energy spectra (from decays of particles at rest) 
and note that there are no z/g produced in this process. The 160 ton detector 
is then used to search for z/g events via inverse beta decay on protons at 
a distance of 30 m from the neutrino source. The experiment detected an 
excess of 87.9 ± 22.4 ± 6.0 events corresponding to an oscillation probability 
of 0.264 ± 0.067 ± 0.45%. (Note that such an appearance experiment affords 
access to very small mixing parameters.) The observed spectrum of events 
is shown in Fig. 10. Other experiments, especially the KARMEN accelerator 
experiment [77] and the Bugey reactor experiment [78], rule out much of the 
allowed region of parameter space but there is a small region remaining at 
90% confidence in the mass range 0.2 < Am^ < 10 eV^, indicating a minimum 
mass of iriy > 0.4 eV (see Fig. 12). It is also significant that the KARMEN 
experiment studied a shorter baseline, which seems to rule out the possibility 
that the Ue are produced at the source. 

Given the well established values of Am^ in Table 1, it is not possible to 
form a third value of Am^ consistent with the LSND data (note that Am^g + 
Am^a + Am|^ = 0). So one would need to either break CPT invariance (allow- 
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Fig. 10. Observed spectrum of positron energy from beam-induced events in the 
LSND experiment [76]. The black dots are the experimental data, the lower his- 
togram is the estimate of various backgrounds, the middle histogram includes the 
estimated contamination from the source, and the top histogram includes the 
neutrino oscillation hypothesis. 

ing niy 7^ w^p), or invoke additional "sterile" neutrinos that do not have the 
normal weak interactions. In addition, the LSND range of Am^ is marginally 
at variance with recent studies of the cosmic microwave background (see be- 
low). Therefore, it is of great importance to attempt to obtain independent 
verification of this result (see J^.l below). 



3.2 Direct Mass Measurements 



Beta decays with low endpoints, in particular tritium {Q = 18.6 keV), have 
been used for a long time in attempts to measure or constrain neutrino mass. 
(^^^Re with Q = 2.5 keV has been also explored recently [79], but the sensi- 
tivity is still only 21.7 eV, a factor of about ten worse than for tritium.) 

There are several difficulties one has to overcome to reach sensitivities to small 
neutrino masses in the study of the electron spectrum of P decay. Obviously, 
the mass sensitivity can be only as good as the energy resolution of the electron 
spectrometer. Also, since by definition the energy spectrum vanishes at the 
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Fig. 11. Results of measurements of rrip^ from tritium /3-decay experiments since 
1990 [20], showing the steady improvement in the achieved precision. 



endpoint, the fraction of events in the interval AE near the Q = EQ—nie kinetic 
energy endpoint decreases rapidly, as {/S.E/QY. Finally, since the decaying 
system is a molecule or at best an atom, one has to take into account the 
possibility that the sudden change of the nuclear charge causes excitations or 
ionization of the electron cloud, i.e., the presence of multiple endpoints. 

Although there are as yet no positive results from direct neutrino mass mea- 
surements, these efforts have made (and continue to make) substantial progress 
(see Fig. 11). Curiously, many past experiments have been apparently plagued 
by systematic effects that tended to mimic a negative miP ^. As the ex- 
periments were improved, it seems that these problems have been largely 
overcome. Two recent tritium (3 decay experiments have reached impressive 
results [80,81] quoting 95% upper limits for the effective mass parameter 
m^pf'^ = \jY.i lUeil'^mf of 2.5 eV and 2.8 eV. However, these values are now at 
the level where atomic and chemical effects on the phase space distributions 
are significant, and it seems that some of these experiments [81] still occa- 
sionally observe anomalous structures near the endpoint. This will likely be a 
substantial challenge for future measurements of this type. 

Constraints from direct mass measurements also exist on m^^ and m^^, but 
these are much larger (see Section 2. 7 for explanation) . Given reasonable as- 
sumptions, it seems that the oscillation results discussed above and the con- 
straints from tritium decay would imply that the mass values are far below 
the direct measurements for z/^ and Ur so we do not dicuss them in detail here. 
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3.3 Double Beta Decay 



Over the last decade, the methodology for double beta decay experiments has 
markedly improved. Larger volumes of high-purity enriched materials are be- 
ing utilized, and careful selection of materials along with deep-underground 
siting have lowered backgrounds and increased sensitivity. The most sensitive 
experiments use the isotopes '^^Ge, i°°Mo, "^Cd, ^^^Te, and ^^^Xe. For ■^^Ge, 
the hfetime hmit has reached an impressive value exceeding 10^^ years. (The 
experimental results are hsted in [20], for the latest review of the field, see 
[12].) The conversion of the observed lifetime limits to effective neutrino mass 
values requires use of calculated nuclear matrix elements, and there is some 
uncertainty associated with them. Nevertheless, the experimental lifetime lim- 
its have been interpreted to yield effective mass limits of typically a few eV 
and in ^^Ge of 0.3-1.0 eV. One recent report [82], analyzing the ^^Ge data 
from the Moscow-Heidelberg experiment, claims to observe a positive signal 
corresponding to the effective mass {mpp) = 0.39^o;28 eV. That report has 
been followed by a lively discussion [83,84,85,86]. If this finding were to be 
confirmed then it would be a major advance in our knowledge of neutrino 
properties, and in particular it would not only prove that neutrinos are Ma- 
jorana particles, but it would also strongly indicate that neutrinos follow a 
degenerate mass pattern, where Am^/m^ <C 1. 

3.4 Cosmological Constraints 

In the early universe, when the temperature was T > 1 MeV, the high den- 
sity of particles allowed weak interactions to occur prolifically leading to a 
substantial density of neutrinos. As the universe cooled to T < 1 MeV, these 
reactions became much slower than the expansion rate and the neutrinos de- 
coupled from the remaining ionized plasma and radiation (photons). Much 
later (~ 100,000 years), the universe cooled enough that atoms formed and 
the radiation decoupled from the matter. The cosmic microwave background 
(CMB) is the further cooled (through expansion) relic of this period, and con- 
tains information on the distribution of matter at that time in the history 
of the universe. The presently observed distribution of matter (through high 
resolution galaxy surveys) and distribution of radiation (CMB) would both be 
affected by the presence of massive neutrinos in the early universe. Although 
the power spectra of the CMB and the density fiuctuations are both sensitive 
to massive neutrinos, a combined analysis of both observablcs is especially ef- 
fective in addressing the existence of massive neutrinos [17,87]. Thus compar- 
ison of the power spectrum of CMB with the observed distribution of galaxies 
can provide information on the sum (over all flavors) of light neutrino masses. 
An analysis of the recent WMAP data [88] yields the result Y^f^vf < 0.7 eV 
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(95% confidence). It is interesting to note that this significantly constrains the 
remaining region of parameter space allowed by LSND. (In addition, interpre- 
tation of the LSND result in light of cosmological constraints requires careful 
consideration of issues related to the behavior (e.g. thermalization) of sterile 
neutrinos in the early universe. [89]) However, it has been argued that the 
obtained limit also relies upon input from Lyman-a forest measurements, and 
that a somewhat less restrictive limit should be quoted [90,91] Other recent 
work [92] also obtains a less restrictive limit without the use of strong prior 
on galaxy bias. 



4 Near-term Future 



The recent discoveries and revolutionary breakthroughs in the study of neu- 
trino properties have motivated a new generation of experimental efforts aimed 
at resolving the remaining issues and establishing new launching points for 
future explorations. These near-term plans and proposals are, for the most 
part, initiatives that advance the themes we have emphasized above. Neutrino 
oscillation studies are planned to address the LSND result, higher precision 
measurements of ^23 and Am^g, and determination of ^13. Direct mass mea- 
surements are aimed at reducing (or discovering) the absolute value of neutrino 
masses and establishing the hierarchical nature of the neutrino mass spectrum. 
And future experiments in double beta decay hope to demonstrate the Majo- 
rana nature of neutrinos and constrain (or discover) the mass scale. Finally, 
the ever increasing precision of cosmological constraints from measurements of 
the cosmic micowave background promise to provide tighter limits on neutrino 
masses. 



4.1 MiniBooNE 



A major near-term priority for the field is to resolve the issue of the validity of 
the results of the LSND experiment. To this end, the Mini-Boone experiment 
[93] has been built at FermiLab with the main goal to explore the same region 
of parameter space with higher sensitivity. This experiment uses a new 0.5-L5 
GeV high-intensity neutrino source and a 800 ton mineral oil-based Cerenkov 
detector located at 500 m from the source to search for the oscillation modes 
and at a later stage — > Pg. The projected sensitivity for 2 years run- 
ning at full beam intensity in each mode is shown in Fig. 12. This experiment 
has the potential to test the validity of the LSND results with high sensitivity 
and to explore the possible role of CPT violation by studying both z/g and 
Pg appearance. If an effect is seen the collaboration plans to mount another 
detector at further distance for additional studies. The initial run with the i/,. 
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Fig. 12. Projected sensitivity [93] of the MiniBooNE experiment compared with 
LSND and other previous experiments. 

beam began in the Fall of 2002 and first results on oscillations are expected 
in 2005. 

4-2 Determination of O'^d 623 

The phenomenon of neutrino oscillation observed by Super-Kamiokande in the 
atmospheric neutrino experiment requires further study to precisely determine 
the mixing parameters ^23- There are two long baseline accelerator 

experiments with prospects for results in the near term: K2K in Japan and 
MINOS in the US. 

The K2K experiment utilizes the 12 GeV proton beam at KEK to produce a z/^ 
beam with average energy 1.3 GeV which is directed at the Super-Kamiokande 
detector a distance of 250 km [94]. Use of GPS receivers enables clean selection 
of events with the proper time relative to the beam spill. Careful monitoring of 
the muons in the tt "'"-decay region ensures proper aiming of the beam over the 
long baseline to the detector. In addition, a combination of detectors located 
300 m downstream of the vr"*" production target is used to measure the fiux and 
energy spectrum of the z/^ beam. During the first data run in 1999-2001, the 
experiment detected 56 fully-contained t'^ events, compared to the expected 
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Fig. 13. Allowed regions from analysis of the first results from the K2K [94] experi- 
ment, representing about 1/2 of the expected total future luminosity. Dashed, solid 
and dot-dashed lines are 68.4%, 90% and 99% C.L. contours, respectively. The best 
fit point is indicated by the star. 

80.ll5;4 events without neutrino oscillations. Thus, so far the K2K experiment 
appears to confirm the oscillation interpetation of the observed anomaly in 
the atmospheric neutrino distribution with high probability. By combining 
the observed energy spectrum and the reduced event rate, the allowed regions 
shown in Fig. 13 are obtained. The experimental plan is to continue running 
and roughly double this data sample. 

The MINOS experiment uses the Main Injector proton beam at Fermilab to 
produce a z/^ beam in the energy range 3-18 GeV directed at a large detector 
located in the Soudan mine at a distance of 735 km [95] . The lowest energy neu- 
trino beam will be optimal for studying the region of Aml2 — 2.0 x 10~^ eV^. 
This experiment also uses a near detector (1 kton) in addition to the far de- 
tector (5.4 kton). The large far detector consists of 486 layers of alternating 
steel plates with finely (4 cm) segmented plastic scintillator planes. The steel 
is magnetized by the return flux of a coil located on the detector axis, which 
enables measurement of muon momentum. With 10 kton-years of exposure at 
full luminosity, the MINOS experiment should easily confirm the z/^ disappear- 
ance observed with atmospheric neutrinos and determine the mass parameter 
Am|2 to about 10%. By searching for appearance the MINOS experiment 
will be sensitive to ^13 in the region sin^ ^13 > 0.15 over the presently allowed 
range of Am32, a factor of ~ 2 lower than the CHOOZ limit. The far detector 
is complete and operating. It is expected that the MINOS experiment will 
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receive the neutrino beam and start operation in 2005. 

There is also a long baseline facility from CERN to the Gran Sasso Laboratory 
(L ~ 730 km). This is a higher energy beam {E^ ~ 17 GeV), suitable for pro- 
duction of r's associated with i/^ — > i>t oscillations. There will be two detectors 
at Gran Sasso for this study: OPERA and ICARUS [96]. The OPERA exper- 
iment will employ photographic emulsion to identify the 'kinks' associated 
with the short (few /xm) r decays. The ICARUS experiment utilizes a liq- 
uid argon time projection chamber to kinematically reconstruct the r leptons. 
Both experiments expect 1-5 events per year to firmly establish this oscilla- 
tion mode, and will begin operation in 2006. Future experiments providing 
precision measurements of the corresponding parameters are highly desirable. 



4-3 Studies of 613, Neutrino Mass Hierarchy, and CP Violation 

The role of mixing between the first and third generation neutrino mass eigen- 
states is pivotal in terms of the phenomenological consequences. The possibil- 
ity of CP violation and implications for leptogenesis scenarios in generating 
the matter-antimatter asymmetry in the universe require mixing between these 
two states. Therefore, establishment of non- vanishing mixing (i.e., 613 7^ 0) is 
of paramount importance and substantial experimental efforts are planned to 
address this issue. 



4.3.1 High-precision Reactor Neutrino Experiments 

When completed, the KamLAND experiment will significantly reduce the al- 
lowed region for Am^2 ^an^ $12 relative to the present results shown in 
Fig. 7. The next major goal for the reactor neutrino program will be to at- 
tempt a measurement of sin^^^i3. As can be seen from Eq. (14) and as dis- 
cussed in section 2 above, such experiments have the potential to determine 
6*13 without ambiguity from CP violation or matter effects. The strategy is 
to capitahze on the success of CHOOZ and KamLAND to build a new ex- 
periment at the 1.5 — 2 km distance now indicated by the Super-Kamiokande 
atmospheric neutrino results. Obtaining the necessary statistical precision on 
this challenging disappearance experiment will require large reactor power 
(> 5 GWth) and large detector size (probably ~ 100 ton). Systematic errors 
must be carefully controlled through reduced background and utilization of a 
two detector scheme (probably in a configuration with a near detector at a 
close ~ 100 — 200 m distance from the reactor(s) in addition to the far de- 
tector at ~ 2 km). Backgrounds must be controlled through sufficient depth 
underground (> 300 m to reduce cosmic ray induced spallation products), 
careful choice of detector materials, and passive shielding (e.g., buffer of water 
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or mineral oil) to reduce events due to radioactive decays in the surround- 
ing rock and other material. Comparison of the rates and observed spectral 
shapes in the two detectors will reduce the sensitivity to reactor source un- 
certainties and provide substantial improvements in the sensitivity to reactor 
antineutrino disappearance in this region. With the anticipated statistical and 
systematic uncertainties held to a total of order 1%, the sensitivity will reach 
about sin^ 26*13 - 0.01 - 0.02 at the optimal value of Am^ ~ 2.0 x lO"-'^ eV^ 
derived from the Super-Kamiokande results. Significant efforts are presently 
underway to optimize the experimental design and select appropriate sites for 
this future experiment. [97,98,99] 

Fig. 14 shows the results of a general study of the potential of such experi- 
ments [97]. With reasonable systematic errors (< 1%) for normalization and 
energy cahbration it is apparent that about 400 GW-ton-years of luminosity 
would provide a determination of sin^ 2^13 with sensitivity better than 0.02. In 
addition, the figure shows interesting behavior at higher luminosity (~ 8000 
GW-ton-years) where the sensitivity improves as l/\/£ due to the high sta- 
tistical precision in determination of the relative spectral shapes at the near 
and far detectors. (The curve with anorm and acai — 00 describes the situation 
where the absolute efficiency and energy calibration of the near and far de- 
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Fig. 14. Projected sensitivity (90% C.L.) of future reactor neutrino experiements 
to ^13 as a function of integrated luminosity [97]. The authors assumed 
Arrig-^ = 3.0 x 10~^ eV^ and studied scenarios with various systematic errors for flux 
normalization (cTnorm) and energy calibration (o"cai)- The JHF-SK line is the pro- 
jected sensitivity of the future JPARC-SK experiment (see J^.3.'2, below), Reactor-I 
scenario is 400 GW-ton-years, and Reactor-II scenario is 8000 GW-ton-years. 
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tectors becomes irrelevant, but the relative efficiencies and energy calibrations 
must be still tightly controlled.) Although the realization of siting KamLAND- 
scale detectors near large power reactor plants would be very challenging, the 
increased sensitivity to sin^ 2612, at such high luminosities would be of great 
interest. 

4.-3.2 Long Baseline Accelerator Experiments 

While the reactor neutrino studies at ~ 2 km have the potential to establish a 
non- vanishing mixing and a numerical value for ^13, further studies related to 
the neutrino mass hierarchy and the role oi CP violation will require mount- 
ing new long basehne accelerator experiments (see Sections 2.4 and 2.5). The 
potential for studying these phenomena is evident in Eq. (32) and illustrated 
in Fig. 15 where one can see the effects of varying the CP-violation param- 
eters and changing the mass hierarchy. There is a great deal of activity in 
this subject at the moment, and a variety of proposed scenarios are under 
discussion. 

The main features of the next generation long baseline accelerator experiments 
include: 

1. ) abihty to detect the i/^ — > fe process as the major goal, 

2. ) L/Ei, ~ 500 km/GeV to optimize the oscillation probability for Am^g ~ 

3 X 10-3 eV^, 

3. ) propagation of the v beam through the earth, resulting in sensitivity to 

matter effects, 

4. ) possibility to vary L/Ej, by adjusting the focusing horn, target position, 

and/or detector location, 

5. ) possibility to switch to 

A significant new concept, the "off-axis" neutrino beam [100] (combined with 
a substantial increase in the neutrino Atix and often referred to as a "super- 
beam"), appears to be a very attractive option for these experiments. By 
positioning the detector off the symmetry axis (at a so-called "magic" angle), 
the neutrino energy becomes stationary with respect to the pion energy and an 
essentially monoenergetic neutrino beam is obtained. Although generally some 
loss of flux results from the off-axis geometry, the advantages include a tune- 
able narrow-band beam with significant suppression of the higher energy tail 
and very low contamination. Thus one can selectively scan a range of Am^ 
to measure the dependence of the oscillation probability with reduced back- 
ground rates (both due to the low Vf, flux and the suppressed 7r° production 
from higher energy neutral current events). This method has been adopted for 
both the JPARC-SK experiment and the Fermilab NUMI proposal. 

The JPARC facility [101] will be a high luminosity (0.77 MW beam power) 
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Fig. 15. Oscillation probabilities for neutrino and antineutrino oscillations [103] 
for the NUMI experiment parameters [102]. The different loops correspond to var- 
ious values of sin^ and normal and inverted hierarchy. The variation of the 
CP-violating phase 6 traces the loops. 



50 GeV proton synchrotron with neutrino production facility aimed about 2° 
off-axis from the Super-Kamiokande detector {L ^ 295 km). The resulting low- 
energy (~ 0.7 GeV) neutrino beam at the Super-K site will enable sensitivity 
to sin^ 26'i3 ~ 6 X 10~^ after about 5 years running. It is expected that the 
experiment would start in 2008, and upgrades involving siting a IMt detector 
and increased beam intensity are envisioned for the future, with potential 
sensitivity to sin^ 2^13 < 1.5 x 10"'^ and CP-violating phase 6 ~ ±20°. 

The "Off-axis NUMI" proposal [102] is to site a new detector either near 
Soudan at L ~ 715 km or at a further site in Canada at L ~ 950 km, at an 
angle of about 14 mrad with respect to the beam axis. The neutrino beam 
energy will be about 2 GeV, and with a 20 kt detector (still under develop- 
ment) the sensitivity would be comparable to the JPARC-SK experiment, but 
the higher energy of this neutrino beam improves the sensitivity to the sign 
of Am23 (i-^-' normal vs. inverted hierarchy). 

There is also discussion of a plan [104] to upgrade the AGS at Brookhaven 
to higher beam power and construct a neutrino beam aimed at the proposed 
new underground laboratory (NUSL) in Lead, South Dakota or some other 
underground location at a similar distance. This would afford the opportunity 
to perform measurements with a large ( 0.5 Mt) detector at a distance of 
L ~ 2500 km, as well as with a closer off-axis detector (L ~ 400 km). Such 
measurements could convincingly demonstrate the oscillatory behavior of the 
flavor transformations with high statistics. 
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Fig. 16. Conceptual design for the KATRIN experiment for higher-precision mea- 
surements of the neutrino mass using the endpoint of the tritium /3-decay spectrum. 

4-4 Future Direct Mass Measurements 



A new experiment to study the tritium /?-endpoint spectrum, KATRIN, is un- 
der development [105]. In order to improve the sensitivity to neutrino masses 
by an order of magnitude it is necessary to both reduce the energy resolution 
of the spectrometer (to ~ 1 eV) and increase the tritium source strength (to 
improve statistical precision). The basic strategy for achieving these goals is 
to scale up the previous spectrometer design used in the successful Mainz and 
Troitsk experiments to a larger physical size, as shown in Fig. 16. The larger 
size enables a higher ratio of magnetic fields -Bmax/-Bmin to improve the energy 
resolution and a larger source acceptance to increase the luminosity. The pre- 
vious spectrometers were 1-1.5 m in diameter and the new KATRIN design 
is 7 m diameter. Both a windowless gaseous tritium source and a condensed 
source will be utilized to enable studies of potential systematic effects asso- 
ciated with the different source methods. The goal for the tritium source is 
a column density of 5 x 10^^ molecules/cm^ and a maximum accepted take- 
off angle of ^max = 51°. These parameters will enable sensitivity to neutrino 
masses in the sub-eV range, hopefully down to ~ 0.35 eV. 



4.5 Plans for Double Beta Decay 



In contrast to the future direct neutrino mass measurements described above 
in 4-4: the next generation of double beta decay searches is poised to ad- 
dress the mass scale indicated by the atmospheric z/ measurements of Super- 
Kamiokande (Am|2 — 2.0 x 10^'^ eV^) and perhaps even approach the scale 
of the solar u solution (Am^2 — 1 x 10^^ eV^). For the case of normal hierar- 
chy (mi ~ m2 -C rns), there is a rather firm prediction for the effective mass 
relevant to decay: 

~ lO'^eV^ . (47) 



m 



Max 



Sin ^2^7712 I 
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This mass scale of less than 3 meV is difficult to address in the near future, 
but if nature actually chooses the inverted hierarchy (mi ~ m2 ^ ^3), then 
the 13 [3 prediction becomes 

mlp ~ If/eimiT + ~ ~ 10"^^ (48) 

corresponding to an effective mass scale of about 30 meV. For the degenerate 
neutrino mass pattern (mi ~ 1112 ~ ms ^ -y/Am|^) the effective mass is larger 
than ~ 50 meV, constrained from above by the mass limit from the tritium 
beta decay. 

The relation between the effective mass in (3(3 decay and the mass of the 
lightest neutrino, evaluated for the solar v solution (see Table 1), is shown in 
Fig. 17. 

Present estimates of the nuclear matrix elements [12] involved in the (3(3 pro- 
cess indicate that, with of order several tons of enriched material, experiments 
could reach this interesting mass range. 




0.001 0.010 0.100 1.000 

minimum neutrino mass (eV) 



Fig. 17. Dependence of the effective Majorana mass {mpp) derived from the rate 
of neutrinoless double beta decay {l/T^^ ~ {f^pp)'^) on the absolute mass of the 
lightest neutrino. The stripes region indicates the range related to the unknown 
Majorana phases, while the cross hatched region is covered if one a errors on the 
oscillation parameters are also included. The arrows indicate the three possible 
neutrino mass patterns or "hierarchies" . 
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Table 2 

Proposed or suggested future Oi^PP experiments separated into two groups based 
on the magnitude of the proposed isotope mass". 









Sensitivity to 


range of 


Experiment 


Source 


Detector Description 


T'^^2 (y) 


{mpp) (meV) 


COBRA[106] 


1 '-IMm 


10 kg CdTe semicond. 


1 X 10^^ 


700-2400 


DCBA[107] 


1 "^n-NT 1 

^^"Nd 


20 kg ^"""Nd layers 
between tracking chambers 


2 X 10^^ 


35-50 


NEMO 3[108] 


iooMo 


10 kg of 0ul3l3 isotopes 


4 X 10^4 


270-1000 




(7 kg Mo) with tracking 






CAMEO [109] 


^^'^Cd 


1 t CdW04 crystals 
in liq. scint. 


> 10^*^ 


> (70-220) 


CANDLES [110] 


^«Ca 


several tons of CaF2 
crystals in liq. scint. 


1 X 10^6 


160-300 


CU0RE[111] 


1 Qn m 

ISOrpg 


750 kg Te02 bolom. 


2 X 10^6 


50-170 


EX0[112] 


i36xe 


1 t «"^Xe TPC 


8 X 10^6 


50-120 


GEM[113] 


^6Ge 


1 t '^"'■Ge diodes 
in liq. nitrogen 


7 X 10^7 


15-50 


GENIUS [114] 


^6Ge 


1 t 86% ^"''Ge diodes 
in liq. nitrogen 


1 X 1028 


13-42 


GS0[115,116] 




2 t Gd2Si05:Ce crystal 
scint. in liq. scint. 


2 X lO^'' 


65 


Majorana[117] 


^^Ge 


0.5 t 86% segmented 


3 X 10^^ 


24-77 




^'""Ge diodes 






M00N[118] 


iooMo 


34 t ""*Mo sheets 
between plastic scint. 


1 X 10^7 


17-60 


Xe[119] 


i36xe 


1.56 t of '^"''Xe 
in liq. scint. 


5 X 10^6 


60-150 


XMASS[120] 


i36xe 


10 t of liq. Xe 


3 X 10^6 


80-200 



" Adopted from [12]. The sensitivities are those estimated by the collaborators 
but scaled for 5 years of data taking. These anticipated limits should be used with 
caution since they are based on assumptions about backgrounds for experiments 
that do not yet exist. Since some proposals are more conservative than others in 
their background estimates, one should refrain from using this table to contrast 
the experiments. The range of the effective masses {mpp) reflects the range of the 
calculated nuclear matrix elements (see Table 2 of Ref. [12]). Again, caution should 
be used since for some nuclei, in particular the deformed ^^'^Nd and ^^''Gd, only few 
calculations exist and the approximations are even more severe than in the other 
cases. 

There are many proposed experiments to address this issue in the near future. 
Most are still in the development stage, and of course the issue of backgrounds 
is critical in every case. Since the source mass is ~ 100 times larger than in the 
present experiments, the background per unit mass must be correspondingly 
smaller. Different proposals approach this issue differently, and in all cases 
substantial R&D is required. Nevertheless, it appears that several experiments 
will be mounted during the next decade with the goal of studying Oi//3/3 with 
sensitivity below 30 meV. The Table 2 lists the proposed experiments we are 
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aware of, together with the claimed sensitivity to the Ou/SP halflife, and its 
crude translation into the sensitivity to the Majorana mass. 

Among the listed experiments, a few relatively small scale searches are ac- 
tually running (COBRA, NEM03). Others, particularly the 'ton size', will 
be mounted in stages with prototypes of 50-200 kg approved and funded at 
the present time for several of them (one, CUORICINO, the prototype of 
CUORE, with 40 kg of natural Te, is already running in Gran Sasso). Clearly, 
the final decision as to which of these many ideas will be realized depends on 
the outcome of these prototype installations. Nevertheless, since they are still 
~ 10 times larger than the present experiments, we can expect substantial 
improvement in sensitivity relatively soon. 

4-6 Cosmological Input 

In the near-term future, substantial additional information from cosmological 
observations will become available that will further constrain the sum of light 
neutrino masses X^/^^i^/ compared to the limits quoted in 3.4- Already in 
progress is the Sloan Digital Sky Survey (SDSS) [121] that will map over 
10® galaxy redshifts and provide relative sensitivity in the relevant region of 
the power spectrum of ~ 1%. When combined with the WMAP CMB data, 
this should provide sensitivity to Yjf^uf ~ 0.3 eV [87]. Further in the future, 
higher precision measurements of the CMB anisotropies are expected from the 
PLANCK mission (now expected to be launched in 2007) [122]. An estimate 
of the combined sensitivity {la) of the expected PLANCK data and the SDSS 
observations to neutrino masses gives the result J2f^uf ~ 0.06 eV [123]. 



5 Longer-term Outlook 

It is difficult to envision what direction the study of neutrino mass and oscil- 
lations will take in the longer-term. Thus, we can only offer educated guesses. 

Some of the future research will, obviously, depend on the results of the near- 
term future described in the preceding section. Among the results that might 
force a revision of the experimental program are the MiniBooNE tests of the 
LSND evidence for the — > Pg oscillations corresponding to Srn^ > 0.2 eV^, 
and the attempts to determine the mixing angle ^13. 

If MiniBooNE confirms the LSND observation, many more experiments will 
be needed because the nature of the neutrino mixing matrix would be much 
richer (more than just three active neutrinos or CPT invariance violation) 
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than envisioned so far and summarized in Table 1. Also, if it turns out that 
the mixing angle ^13 is relatively large, perhaps close to the present upper 
limit, the study of the possible CP invariance violation in the neutrino sector 
will become considerably easier and will proceed faster than if ^13 is very small. 

In addition, neutrino astrophysics will undoubtedly advance. Further studies 
of the solar neutrino spectrum will be conducted, and in particular the low 
energy part, the flux of the ^Be and pp neutrinos, will be likely determined 
in live counting experiments. These measurements will be valuable not only 
for the further refinement of the determination of the oscillation parameters 
Amf2 and ^12, but also for better understanding of solar energy production. 

Moreover, with several experiments on-line, we can hope that in the near 
future the next galactic core collapse supernova will be observed by a variety 
of neutrino detectors. If that happens, the observation of SN1987A [124] which 
launched this field of neutrino astrophysics, will be exceeded many times, and 
all components of the supernova neutrino spectrum will be observed. Again, 
such observations would not only advance our knowledge of neutrino physics, 
but will offer a unique insight in the physics of the stellar core collapse. 

We also anticipate that during the next decade other applications of neutrino 
physics will become reality. For example, the observation of 'geoneutrinos', 
i.e. the z/g emitted by the decay of radioactive U and Th series in the Earth 
crust [75,125,126] would enable determination of the corresponding radiogenic 
heat generation and provide an important contribution to geophysics. Another 
application is the measurement of the relic supernova flux: the diffuse flux 
resulting from past supernovae in the surrounding space, roughly up to redshift 
z = 1. By measuring this flux one can determine the average supernova rate 
over a substantial part of the visible universe (for an estimate of the rate see 
e.g. [127]). 

At present, neutrino telescopes have not yet seen the TeV and higher energy 
neutrinos which likely accompany the production of cosmic rays with such 
energies. This might change in near future when the relatively small under- 
ground detectors are supplemented by much larger ones in open water or ice. 
Observation of the direction of high energy neutrinos would considerably ad- 
vance the study of the origin of high energy cosmic rays [128]. It will also make 
it possible to search for neutrinos from annihilation of the so far hypotheti- 
cal dark matter particles (WIMPs). It is beyond the scope of this review to 
describe the projects being built or proposed (see [129] for a recent review). 

Naturally, the whole field of particle physics will advance at the same time. 
With the launch of LHC one can imagine that the existence of the Higgs boson 
will be experimentally confirmed, and its properties will be, at least crudely, 
determined. One can also hope that the quest for the understanding of the 
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two main paradigm of present day physics, the nature of the 'dark matter' 
and 'dark energy', responsible for most of the energy density of the Universe, 
will be advanced. 

These more general areas of particle physics could be, in fact, intimately re- 
lated to the quest for neutrino mass and oscillations. Taking the see-saw for- 
mula AiL = M.D^ / -M-H, and using for the neutrino mass scale ~10 meV 
and for M.£, the electroweak symmetry breaking mass scale ~ 100 GeV, we 
arrive dX M.^ ^ 10^^ GeV, i.e. close to the GUT unification mass scale. If that 
relation could be firmed up, we will arrive at another determination of that 
important mass scale. The relation of the neutrino mass scale to the 'dark 
energy' is based on the, perhaps accidental, numerical coincidence. Remem- 
bering that the energy density of the 'dark energy' is about 0.7f2c ~ 3.5 x 10^ 
eV/cm^ and rewriting it in the 'natural units' as e^/ {hcY we arrive at the dark 
energy mass scale e ~ 2 meV, close to the neutrino mass scale, and unlike any 
other mass scale in particle physics. (For recent attempts to relate these two 
mass scales see e.g., [130].) 

Future large scale projects to generate powerful new beams of neutrinos are 
already envisioned. To explore the possibility oi CP violation in the lepton 
sector one would like to have a well-understood and collimated neutrino beam 
of a well defined flavor and energy spectrum. Such a beam could be aimed at a 
distant large detector. Recently it has been pointed out that technology exists 
to construct pure and Ve beams of the required properties. Accelerating 
radioactive ions, ^He {Eq = 3.5 MeV, T1/2 = 0.8s, produces z/g) or ^^Ne (£'0 = 
4.45 MeV, T1/2 = 0.8s, produces Ue) to 7 ^ 100, would produce beams of 
precisely known energy profile extending to 2'jEq, and collimated to I/7. The 
oscillation signature in the far away detector would be the appearance of /i"*" 
or /i~ [131]. It is expected that 0(10^^) decays per year can be achieved. 

The beams of z/^ and as well as z/e and neutrinos could be obtained in 
a neutrino factory where accelerated muons are stored in a ring with long 
straight sections. Such a beam will produce neutrino beams with equal mix- 
tures of z/^ and if /x"*" are stored, and with if /i~ are stored. It is 
expected that neutrino factories could provide O(10^°) useful muon decays 
per year. Neutrino factories would thus provide beams with small systematic 
uncertainties in the beam flux and spectrum. CP violation, and determina- 
tion of sin^ 29i3 even if it is as small as 10~^ — 10~^ can be achieved with such 
beams and baselines of several thousand km [132]. 

One indication of the importance of the study of neutrinos is the recent report 
"Facilities for the Future Science, A Twenty- Year Outlook" (see "http://www.se. doe. gov/ I. 
It lists two projects, among the 28 listed in all fields supported by the US De- 
partment of Energy, relevant to the present topic, the study of double beta 
decay in underground detectors treated among the mid-term projects and Su- 



47 



per Neutrino Beam treated as a far-term project. 

Clearly there are high expectations in the community of particle physicists 
that experiments with neutrinos will continue to provide new insights and 
substantial progress in the coming decades. Construction and operation of 
the required facilities will be challenging, but the potential rewards are great, 
and a large and talented community of physicists is poised to embark on this 
unique and interesting journey. 
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